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Abstract 

The problem of heterogeneous nucleation of second-phase in alloys in the vicinity of elastic defects 
is considered. The defect can be a dislocation line or a crack tip residing in a crystalline solid. We 
use the Ginzburg-Landau equation to describe the spatiotemporal evolution of the order parameter 
in the environs of the defect. The model accounts for the elasticity of the solid and the interaction 
of order parameter field with the elastic field of the defect. A finite volume numerical method is 
used to solve the governing partial differential equation for the order parameter. We examine the 
nature of the phase transition in the vicinity of the defects. 
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I. INTRODUCTION 



The presence of elastic defects such as dislocations and cracks may induce nucleation 
of a second phase in many alloys . For example, the formation of brittle hydrides in 
titanium and zirconium alloys (TiH-,,, ZrH x ) is of special interest for aerospace and nuclear 
industries, since they may cause embrittlement of these alloys used in various equipment jsj. 
Hydride formation is commonly accompanied by a preferred orientation of the precipitates 
(platelets) due to the texture of the polycrystalline material and/or the presence of external 
stress Q]. Moreover, the crystal structure of the hydride (face-centered cubic for 5-hydride) 
differs from that of the matrix (hexagonal close-packed for a-Ti). In the present study, we 
use a phase field approach to analyze the nature of the structural part of the new phase 
formation near elastic defects. That is, the effect of composition is not included in our 
analysis. The model used here rests on the Ginzburg-Landau theory of phase transition in 
which a scalar non-conserved order parameter characterizes the presence or the absence of 
the nucleus. The interaction between the order parameter and the deformation field is also 
taken into account [5, 6|. A more general set-up was presented in 



II. MODEL DESCRIPTION 



We consider two types of elastic defects in a crystalline material; an edge dislocation and a 
semi-infinite crack. The considered phase transformation is the nucleation of a second phase 
in an elastic material. A single structural order parameter that accounts for the symmetry 
of structure is used to characterize the phases. It is supposed to be a scalar field (Ising 
model) r](r,t) that is a function of space r and time t. Hence, rj — corresponds to solid 
solution and rj to a nucleus. The total free energy of the system is written [3] 

T — T s t + T e l + Tint-, (1) 

where T s t is the structural free energy, T e i the elastic strain energy, and Tint is the interaction 
energy between the structural order parameter and the strain field. The structural free 
energy is 

= J [|(V77) 2 + V(77)]dr, (2) 
2 



where the space integral is within the volume of the system dr = pdpdOdz. Here g(Vr]) 2 
accounts for the spatial dependence of the order parameter, g is a positive constant, and the 
second term in the integrand is the Landau potential 8] 



(3) 



where r$ is taken to be a linear function of temperature T, e.g. tq = ao(T — T c ), ao is a 
positive constant and T c the transition temperature in the absence of elastic coupling. The 
coefficients uq and vq are considered to be temperature independent. The elastic free energy 
is 

2n 



T A = J [f(V-u) 2 + M^(n..-^V- U 



dr, 



(4) 



where K and M are the bulk and shear modulus, respectively, = (Vjttj + Viii 3 -)/2 is the 
strain tensor with V< = d/dxi, d the space dimensionality, and i,j stand for x, y, z in d = 3 
(x, y in d = 2). Finally, the interaction energy is 



int 



T] 2 V • udr, 



(5) 



where rfV ■ u describes the interaction between the order parameter and the deformation. 
The strength of this interaction is denoted by k and is taken to be a constant. 

The temporal evolution of the spatial order parameter is determined by solving the time- 
dependent Ginzburg-Landau equation for a non-conserved field, cf. {9], 
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drj 
dt 



Sr) 



(6) 



where L a is a kinetic coefficient that characterizes the interface boundary mobility. Defects 
in a crystalline material, such as dislocations and cracks, render internal strains which change 
the equilibrium condition in the solid. If f (r) denotes the variation of the strain field due to 



the defect, the equilibrium condition that includes the force field is generated by 

MV 2 u + (A - M)VV ■ u + nVrf = Mf (r), 



10| 



(7) 



where A = K + 2M(1 — 1/d). Equation (J7|) is then used to eliminate the elastic field from 
the expression for the total free energy, which now can be expressed as 



dr. 



(8) 



The last three terms in Eq. (j8]) correspond to the Landau potential energy in Eq. (3) but 
with modified coefficients: 



r\ = r — kAcos9/p, for an edge dislocation, 
ri = r — kB cos(6 l /2)/p 1 / 2 , for a crack, 
ui = uq — 2k 2 /A, 



(9) 
(10) 

(11) 



with A = (26/vr)M/A and B = 4fC T (l - 2u)(l + v)j (2n) 1 / 2 E, where b is the magnitude of 
the Burgers vector, K\ the mode I stress intensity factor, v Poisson's ratio, and E Young's 
modulus of the material. For a defect free crystal r x = r , and for a rigid crystal U\ = Uq. 
Figur^H shows the geometry of the two defects. 



(a) Edge dislocation (b) Semi-infinite crack 

FIG. 1. Geometry of the defects; x = pcosO and y = psinO. 



Inserting Eq. (|H]) into Eq. ((6]), gives us the governing equation for the space-time variation 
of the order parameter 

1 dr] 



L a dt 



(12) 



At equilibrium, drj/dt = 0, and the Landau potential, i.e. Eq. ([3]), tells us the kind of 
phase transition plus for which set of parameters nucleation and growth will take place. 
The Landau potential with the modified coefficients r\ and U\, and v o > for two sets of 
parameters are shown in Fig. [21 one with u\ > and the other with u\ < 0. If u\ > 
and r\ > 0, no nucleation of the second phase will occur, however, r\ = embodies an 
onset of nucleation, and for r% < the nuclei will continue to grow since only the second 



phase is stable, see Fig. 2(a) In the case of ui < (Fig. 2(b)), a metastable second 



phase may emerge if r± < uf/Avo since the potential then has two local minima at 77 7^ 0. 
For r\ = 3uf/lQv , both phases are equally stable and so may coexist. For even a smaller 
positive value of r± the solid solution is metastable and the second phase stable. When 
r\ = only the second phase will exist. 




(a) ui > (b) mi < 

FIG. 2. The Landau potential V{rf) = \rirj 2 + |ui^ 4 + gfo?? 6 with vq > 0. 



Since r\ is a function of both temperature and spatial coordinates, the location where 
nucleation and further evolution of the second phase will occur, relative to that of the defect, 
will vary with the choice of U\. In Fig. [31 the boundaries corresponding to triple minima, 
i.e. r\ = 3u1/16 when u\ < 0, for the dislocation and the semi-infinite crack are shown, 
respectively. In the figure, normalized coefficients U\ and R\ corresponding to u\ and r\ are 



used, which are defined below. In Figs. 3(a) and 3(c) with R > 0, there are areas close to 



defect where Ri < 0, and thus a second phase is expected to nucleate there for all values of 



III. NUMERICAL METHOD 



For the numerical simulations, we have used the open-source program FiPy ll| which 
is an object oriented, partial differential equation solver based on a standard finite volume 
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(c) Rq > (d) Rq > 

FIG. 3. Characteristic contours in the presence of a dislocation (a)-(b), and a crack (c)-(d); where 
R Y = 3E7?/16 for C/i < 0. 



approach. A square 2D-mesh consisting of 200 x 200 equally sized elements is used for the 
computations, where each element has a side length A/ = po/20, with po being a charac- 
teristic length that is introduced below. The simulations are performed using a constant 
time increment At = 0.9AZ 2 / (2g). Periodic boundary conditions are applied, and the initial 
value of the order parameter is taken to be a small random number between 0.005 and 0.01. 
The dislocation is located at the origin, (x, y) = (0, 0), with the slip direction along the line 



x = 0, see Fig. 1(a) The crack lies along the negative part of the x-axis with the tip at the 
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origin, see Fig. 1 (b) The system is considered to be large enough to model a single disloca- 
tion and a long single crack, respectively, without being affected by the periodic boundary 
conditions. For the computations, Eq. (IT21) is rewritten as 

|Z = V 2 r/ - (R lV + U lV 3 + V 5 ) , (13) 

where t = L a vot, x = (vo/gY^x and y = (vo/gY^y. The coefficients are now i?i = 
i?o(l — Po cos6/p) for the dislocation with R = Tq/vq and p — kA/tq] whereas R\ = Rq[1 — 
cos(6'/2)/(p/po) 1 ^ 2 ], Po = (K-B/r Q ) 2 for the crack, and U\ = ui/vq. Thus a characteristic 
length, po, is introduced. In this parametric study, Ro and XJ\ are varied and p is set 
equal to unity. The time steps are taken to be small enough to capture the evolution for 
all the different combinations of Rq and U\ that are studied. Rq is a linear function of 
the temperature difference, as defined above, and thus quenching the system to below the 
transition temperature means R Q < 0; and for T > T c , R > 0. 

IV. RESULTS AND DISCUSSION 

Results from simulations of the evolution of in the vicinity of a dislocation are first given 
for some combinations of Ro and U\. Thereafter, the outcome from the calculations with a 
crack in the system is alluded. The details for the latter case will be presented elsewhere. 

A. Dislocation: Rq = — 1 and U\ = 

The case Ro — 1 and U\ = corresponds to a quick decrease of temperature, i.e. quenching 
below the transition temperature. In Fig. HJ the evolution of t] is illustrated by surface plots 
representing four different times. One can clearly see that a top is growing below the 
dislocation (x < 0) with its peak close to the dislocation. The evolution of t] is also shown 



in Fig 5(a) as profiles of rj (y = 0) at different times (t = 50 — > 650At in steps of 50 At). The 
top first grows until the peak reaches a stable value (here ~ 1.5), while in the surroundings 
7] pa . Thereafter, the top broadens and 77 increases also for x > 0. However, close to 
the dislocation, the increase is held back. This process will continue until rj ~ 1 covers the 
whole material, leaving a top and a valley on each side of the dislocation, respectively; see 



the contour plot in Fig. 5(b) The dashed circle represents i?i = 0, and inside it R 1 > 



which, according to the Landau potential, indicates that no nucleation should take place 



there, cf. Fig. 2(b) Since the magnitude of R\ goes to infinity as while approaching 



the origin, the peak and the valley are expected to evolve. Far from the dislocation, Ri goes 
to i?o, i-e. a constant negative value, and therefore nucleation is expected. 
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FIG. 4. Spatial distribution of 7] at various times t = 100, 150, 200, 250At, for Rq 
AT < 0), and U\ = 0. Only every fifth node in the mesh is used for the illustration. 
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FIG. 5. R = -1 and U t = 0: (a) Evolution of r](y = 0) for t = 50 -)• 650 At. The dotted vertical 
lines indicate p = ±po! (b) Contours of r/ at t = 650 At. At the dashed circle the sign of R\ changes 
and is positive inside the circle. 



B. Dislocation: Rq = — 1 and Ui = — 1 

Figure 6a shows the evolution of r](y = 0) whereupon the same trend as for XJ\ — is 
observed. First a top is growing near the dislocation, where an increase of r\ arises from this 
top to finally form a plateau. However, for Rq = —1 and U\ = —1, the evolution is faster 
than in the foregoing case, and both the peak and plateau values are higher, i.e. r] = 1.7 and 
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1.3, respectively. It should be noted that the time increments are equal for all the studied 
cases. 
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(b) E7i = l 

FIG. 6. Evolution of r/(y = 0) for t = 50 — s> 650Ai, upwardly, with f?o = —1. The dotted vertical 
lines indicate p = ±po. 



C. Dislocation: i?o = — 1 and U\ = 1 

In the case i?o = — 1 and U\ = 1, the development of 77 has the same character as the 
foregoing cases, with the difference that the maximum value of 77 and the plateau level are 



smaller and the evolution is much slower. In Fig. 6(b) only profiles up to t = 650 At 



are shown, however, for the evolving plateau, 77 reaches around 0.8. It is observed that 
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the plateau values for the three different cases correspond well to the real, positive roots 
rj = (1.00,1.27,0.78) that give the minima of the Landau potential for Ri = Rq = —1. 
These non-trivial roots are found from 



2 -U x + y/U? - 4R Q 



vi = 1 v 1 (14) 



D. Dislocation: Rq = 1 

Now consider the situation where the temperature is above the transition temperature 
(Ro > 0), i.e. no phase transition is expected in a non-stressed and defect free material. 
However, by introducing a dislocation into the system the situation alters. Figure 7 shows 
the evolution of r\ in the case where U\ = (—1,2, —3), repectively. For U\ — — 1 a top emerges 
and evolves until it finds a stable shape, with U\ = —2 a top develops and thereafter a small 
plateau appears that will stop to broaden at t ~ 2000At, and for U\ = —3 the top grows 
and broadens until the plateau covers nearly the whole area except near the dislocation 
where x < 0. In Fig. 3a the contours where both phases have minima of the Landau 
potential are shown for different values of U\, and the more negative is U\, the larger is the 
circular area where the second phase is assumed to evolve. For U\ = —1, the diameter of the 
corresponding circle is approximately 1.2 and for JJ\ = —2, the diameter is 4.0. Dashed lines 
in the figures show the corresponding location on the a;- axis. In the case with U\ = —3, the 
contour surrounds a small circular area below the dislocation (x < 0) wherein no analogous 
phase transformation is expected. 



E. Crack 

Calculations have also been performed for the crack and the results are similar to that 
of the dislocation. In the vicinity of the crack, there is always an area where the second 
phase is stable regardless of the sign of Rq. Detailed results and analysis for the crack will 
be presented elsewhere. 
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FIG. 7. Evolution of rj(y = 0) for t = 50 — > 650Ai, upwardly, with Rq = 1. The dashed curves 
represent the r] profile at t = 400Ai. In (b) and( c) profiles for t = 2000Ai are also added (dash 
dot curves). The vertical dotted indicate x = ipo and the vertical dashed lines show the location 
with triple minima. The dashed curves represent the r] profile at t = 400 At. 



V. CONCLUSIONS 

Numerical simulations of the spatiotemporal evolution of a second phase in the vicinity of 
elastic defects in crystalline solids have been performed using the Ginzburg-Landau equation 
for a single non-conservative structural order parameter. The computations indicate that 
these defects always trigger a nucleation of a second phase. In the very vicinity of the 
dislocation and at the crack tip a distinct top emerges and evolves. In some cases, the 
structural order parameter evolves into a plateau, which either finds an equilibrium shape 
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or spreads out into the whole material. It should also be mentioned that the Landau type 
energy could also be used to roughly estimate the spatial evolution of the order parameter in 
the vicinity of defects. The study does not include the influence of concentration of species 
which is essential for modeling phase transition in an alloy. With a two component order 
parameter field the local orientation of the second phase could also be deduced. Analysis 
comprising the non-conserved order parameter coupled to a conserved concentration field 
obeying a diffusion-like equation will be presented elsewhere. 
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